FT-031 



A simple variational approach for an interacting 
Fermi trapped gas. 
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Quantum states of a two-component Fermi trapped gas are described by introducing an effective 
trap frequency, determined via variational techniques. Closed expressions for the contribution of 
a contact interaction potential to the total energy and the pairing interaction are derived. They 
are valid for both few and large number of particles, given the discrete nature of the formulation, 
and therefore richer than the continuous expressions, which are perfectly matched. Pairing energies 
within a shell are explicitly evaluated and its allowed values at a given energy level delimited. We 
show the importance of the interaction over the trap energy as the number of particles (N) grows and 
the temperature decreases. At zero temperature we find a polynomial dependence of the interaction 
energy on the Fermi energy, whose dominant term at large N corresponds with the mean field 
approximation result. In addition, the role of the strength of an attractive potential on the total 
energy is exhibited. 

PACS numbers: 03.75.Ss, 05.30.Fk, 71.10.Ca 



I. INTRODUCTION 

In the last years, trapped degenerate two-component 
Fermi atomic gases have received great interest both from 
the theoretical 

0, H H II H Q and the experimental 00 
point of view. One of the reasons is that these systems of- 
fer an ideal scenario to study pairing interactions at vari- 
able interaction strength, density and temperature. It 
has been recognized that atoms like 40 K and 6 Li exhibit 
magnetically tunable Fcshbach resonances so that the s- 
wave scattering interaction may be varied from repulsive 
to strongly attractive. The accessible states for a di- 
lute gas of trapped atoms may be modeled in the zeroth- 
order approximation as a system of noninteracting par- 
ticles. The simplest theoretical approach for describing 
trapped Fermi atoms when the interactions are not negli- 
gible, consists in approximating the system wave function 
as properly symmetrized products of one-particle eigen- 
states of an effective harmonic oscillator, whose frequency 
is determined via a variational procedure. The purpose 
of this work is to make an study of the predictions of such 
an approach. In particular, we obtain closed expressions 
of the total interaction energy that reflect its increasing 
relevance as the number of particles in the system is aug- 
mented, as well as closed expressions for pairing energies 
of atoms within a given shell. 



II. FORMALISM 

Let us consider a gas of Fermi atoms of mass M with 
internal states | f) and | |). They are in a harmonic 
trap that for simplicity will be assumed to be of spheri- 
cal symmetry and characterized by a frequency u>o- An 
asymmetric harmonic trap can be treated in a completely 
analogous way. We model the interaction between atoms 



in different internal states by an effective contact poten- 
tial 



X,-T . 



(1) 



with cvo the scattering length. This potential is based 
on scattering theory at low energies. It is inexact for at 
least two reasons: it takes into account neither higher 
partial waves nor the dependence of s-wave scattering on 
the momentum transfer. One consequence of the latter is 
that no parameters beyond the scattering length such as 
the effective range are included in Eq. (|TJ . Thus, this po- 
tential is expected to be valid for low energy transfer. In 
fact it has proven to be useful to describe cold atoms at 
low densities Q and it could also be important for inter- 
mediate densities by introducing an effective scattering 
length [h}. 

The Hamiltonian of the system is taken as 



(2) 



We are interested in obtaining approximate eigenfunc- 
tions of this Hamiltonian using the variational method. 
The validity of such an approach is limited by the relia- 
bility of the model Hamiltonian. As a trial wave function, 
we consider the product of two Slater determinants, one 
for each internal state. These determinants are built up 
by the one particle eigenfunctions of a spherical harmonic 
oscillator of frequency to' and quantum numbers n x ,n y 



and 



The value of lu' should reflect the interaction 



effects and in the limit of zero interaction it becomes 
lu' = luq. Notice, however, that dramatic changes in lu' 
with respect to loq are incompatible with the assumptions 
supporting the use of the contact interaction potential 
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Taking into account the results summarized in the Ap- 
pendix, the expectation value of the Hamiltonian for this 
trial function is found to be 



(H)=fkJ [ 
where a = u>' /u>o, 



e* + £o 



a' 
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(3) 



For a given quantum number n, the integrals I(n, m) 
achieve its maximum value for n = m. In Figure 1, we 
plot the allowed region of the 8 n ^ values, delimited by 
solid line curves. They correspond to its maximum S max 
(top) and minimum <5 mi „ (bottom) values, for < n v < 
200. Systematically 5 m ax corresponds to two quantum 
numbers equal to zero, for instance, n y = 0, n z = 0, 



+ nl + 4)+ Yl ( n 
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+ n y + n z ) (4^j ecom p 0se d j e . 

n z <~ l/5n v . 



Meanwhile S m i n is achieved when n v is 
5., in the form n x ~ (3/5)n„ and n y = 



(5) 
(6) 



III. GROUND STATE. 
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min{n,m} 



I(n, m) 
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As a particular application of these formulae we con- 
sider the ground state wave function corresponding to a 
configuration where all the possible states up to the non 
interacting Fermi energy €f — hu n (M f + 3/2) are occu- 
pied. To characterize the system under these conditions, 
(2(m — k) — 1)!! (2(n - k) - 1)!! (2k w - F^f^ to calculate the summations 



with N( % > the number of particles in a given internal state 
the scattering length in natural units and 



fe=0 



(2(m-fc))!! (2(n -*;))!! 

(7) 

e t , eo and e s can be identified with the trap, zero point 
and interaction energies in units of Hujo, respectively. 
Taking a as a variational parameter, the condition that 
the expectation value of the Hamiltonian becomes an ex- 
trcmum reduces to 



(2kfyV. F M F -n x M F -n x -n y 

EE ^ 1 



=0 n v —0 n z =0 
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n v = 



Mi 
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n z =0 



n x = 



{M F + l) 3 + 3(M F + l) 2 + 2(M F + 



(M F + I) 4 + 2(A4jr + l) 3 - (M F + 
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0. 
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If the scattering length is negative, in order to admit a 
variational solution corresponding to a minimum of the 
energy it is necessary that 



The first sum let us evaluate the number of particles 
for such a configuration in terms of the Fermi energy. 
While the second gives one third of the corresponding 
trap energy in fko units, neglecting the zero point 
energy. Notice that for large M.? the trap energy be- 



comes 



7 5/2 



< - 



4(e t + e ) 
3e s 
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Once the proper value of the variational parameter has 
been determined, the properties of the system can be es- 
timated. In particular, the pairing energy between two 
particles of quantum numbers n\ whose wave func- 
tions are (j) n t (x) and 4> n i (x) is given by 



(V) 



d 3 x</)* n T (x)0„ T {x)4>* n i {x)(p n i (x) 



(10) 



AlTCtr, 



=I(n l x ,nl)I(n l y ,nl)I(n l z ,nl). 



For atoms in a trap with effective frequency u)', an energy 
level e nv = hu>'(n v + 3/2) exhibits a degeneracy (n v + 
2)(n v + l)/2 arising from the condition n x + n y + n z = 
n v . Then the pairing energy between two particles in 
the same shell is directly given by the restriction n\ — 
n x + n y +n\ = n* + n y + n\ = n*. In order to illustrate 
the pairing energy features, let us explore the behavior 
of the universal dimcnsionless quantity 



That is, the discrete summation relation between these 
variables coincides with that obtained in the standard 
procedure that works with an approximate continuum 
spectra. In Table I, the behavior of the dimcnsionless 
quantities iVW — — N^, ef* and e s = e s /g as a func- 
tion of A4f is illustrated. Notice the increasing relevance 
of the interaction term e s compared to the e t . While at 
Mf = 10 it represents only a 13.6%, at = 300 it 
becomes 47% as big as et- In fact, for < A4f ~ 200, 
we have found that e s has a polynomial behavior of the 
form 



= (2tt)- 3 / 2 +M 1 J 2 (c +c 1 Mf- 



c 2 M 



with 

co = 
C3 = 



0.11671733, 
0.03706759, 



Cl 

c 4 



0.18747855, 
0.00411862 



■C 3 M F +C 4 M F ) 

(15) 



c 2 = 0.12399665, 



= I{n xi n x )l{n yi n y )I(n z ,n z ), 



n x +n y +n z wrtlma.relative error of the order of 10 7 for A4 F > 40 and 
(11) of order less or equal to 10~ 4 for Mf < 40. That is, our 
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scales as Mp 5 
M% 



numerical calculations show that the interaction energy 
for large Aip growing faster than e t ~ 
This numerical result can be understood within 
the mean held approximation as follows: 

In general, for a given mixture of Fermi atoms de- 
scribed by the product of two Slater determinants formed 
by one-particle eigenfunctions (pi (fj ) of an effective one 
particle Hamiltonian, h = p 2 /2M + v(r), the expectation 
value of a contact interaction potential is determined by 



; (s(xn - x T )) 



TABLE I: Behavior of the dimensionless quantities: number 
of particles N, trap energy e t and interaction energy e a — e s /g 
as a function of Mf- 



<fxp\x)p l {x) (17) 



where p^'(x) represents the spatial density of atoms of 
z-type. In our problem, the one particle Hamiltonian 
is that of a spherical harmonic oscillator of frequency 
Lo' . Then, at zero temperature and for a large number 
of atoms, the Thomas-Fermi approximation leads to the 
expression 0,0] 



P {i) (x)t=o 



ATW 8 r 



R F n 2 
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So that, 
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(18) 
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V / 2tt 3 2835 
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Eq. (|15|l is consistent with the later result and, in fact, 
generalizes this formula for small M. f ■ 

Another useful result to exploit the variational ap- 
proach, for large iV (T) = iV (i) = iV w , is that the varia- 
tional equation © can be approximated by 



0Al9372gR F a 5 - (a 2 - I) 2 = 0, 



(21) 



where R F = R F /y/fr/Muj' = (48ATW) 1 / 6 is the charac- 
teristic size of the trapped degenerate Fermi gas in the 
natural units of the trap. In Table II, the results of the 
variational calculation are exemplified for M. F = 150, 
corresponding to about half a million atoms of a given 
internal state. An attractive interaction is considered and 
the corresponding coupling constant g is varied, leading 
to optimal values of a via the variational equation. No- 
tice that the information about the trap parameter (ujq) 
and the atom characteristics (M and ao) are contained 
in g. Therefore the results are not restricted to neither 
a specific atom nor trap. In this table, the exact expec- 
tation values of the Hamiltonian in Hu>q units divided by 
the total number of particles N = N' + N* are shown 
for lo' — ljq ( (E)u'=u ) and for an optimized choosing of 
the variational parameter a ( (E) m i n ). As predicted by 
Eq. ijjjj no minima can be found for \g\ > 0.9. In this case 
the expectation value of the energy exhibits a maximum 



Mf 




4 l) 







00 


91 n 


9^ A 


1U 




91 A ^ 


zyo.u 


on 


1 771 
1 / / 1 


ZDODO 


/ino o 
401Z.O 


ou 


0400 


1 997^0 


Z4oyo.y 


40 


12341 


370230 


82991.0 


50 


23426 


878475 


217041.8 
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479024.0 
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3265290 


938926.7 
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91881 


5512860 


1685871.8 


90 


129766 


8759205 


2829605.9 


100 


176851 


13263825 


4501889.6 


150 


585276 


65843550 


27103302.6 


300 


1373701 


206055150 


97458248.4 



(19) 



at a given value of a < 1 and decreases monotonically for 
a > 1. That is, the gas would behave as trapped in a har- 
monic potential of arbitrarily large frequency. Therefore, 
it would collapse. However, this prediction is beyond the 
validity region of the model. 

A mechanical instability of a strongly interacting Fermi 
gas arises when attractive interactions overcome repul- 
sive effects such as the Fermi pressure 0- In actual ex- 
perimental realizations [10|, there are mechanisms not 
represented by the simple model Eq.Q that may stabi- 
lize the system. For instance, repulsive contributions to 
the interaction arising from higher order partial waves 
|ll| , or within the s-wave scattering scheme, effective po- 
tentials resulting from adding the ladder diagrams in a 
many-body perturbative approach 0. Assuming densi- 
ties that guarantee an interatomic separation larger than 
the range of the potential but much shorter than the 
scattering length, Heiselberg found that an effective s- 
scattering amplitude arises, which always favors the sta- 
bility of a two component Fermi gas. 

A mixture of spin-polarized atomic 6 Li in two possible 
hyperfine states constitutes a specially interesting sys- 
tem 0, H, El- It exhibits magnetic tunable Feshbach 
resonances that enable the variation of the scattering 
length from strongly repulsive to strongly attractive. In 
fact, an anomalously large and negative scattering length 
ao = —2160 ao has been already measured The 
coupling constant corresponding to this value of ao 5 for a 
typical trap of mean frequency ujq ~ 3kHz is g — —0.725. 
Then the effective frequency for half a million atoms is 
(J = 1.48^0- It differs significantly from the trap fre- 
quency; the model potential and our variational approach 
are at the border of its applicability. 
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TABLE II: Average energy per atom for u)' = luq, 
(E) w i =ulQ J 'N^ and for an optimized variational parameter a, 
as a function of the strength of an attractive 

potential, g . 



g 




a 




-0.05 


111.34 


1.02 


111.33 


-0.10 


110.18 


1.03 


110.13 


-0.20 


107.87 


1.07 


107.63 


-0.30 


105.55 


1.11 


104.99 


-0.40 


103.24 


1.17 


103.24 


-0.50 


100.92 


1.23 


99.12 


-0.60 


98.61 


1.32 


95.79 


-0.7251 


95.71 


1.48 


91.03 


-0.80 


93.98 


1.65 


87.64 


-0.85 


92.82 


1.90 


84.93 



IV. INTERACTION ENERGY FOR TV 0. 

A detailed description of the behavior of the Fermi gas 
at finite temperature including interactions, can be given 
by evaluating the variational wave functions for a set of 
accessible quantum states and performing the proper sta- 
tistical average to infer properties like average energies or 
densities. In this section we shall make an approximate 
calculation where the starting point is Eq. (|16fl . This 
approach has not a rigorous theoretical basis but should 
give us a reasonable estimate of the energy shift at finite 
temperatures. According to the Thomas-Fermi approxi- 
mation, valid for a large number of particles, the atomic 
density for a trapped Fermi gas, at a given temperature 
T, is @ 



pV>%T) = 



1 

2^ 



d 3 k 



3 /3[R|/c| 2 /2M+l/2Mw' 2 r 2 ]-2 i 1 



(22) 



where (3 = 1/ksT. Defining the function 



1(0 = J dn J drf J dCW 2 C 2 [e ,,2+c2 - ? +irV' 2+C2 ^- 

(23) 

we can estimate the interaction energy at such tempera- 
ture by generalizing equation (|16fl to consider tempera- 
ture dependent densities: 

(V(N,T)) = f d 3 xp r (x,T) Pl (x,T) (24) 



M 

47ra V / M2 9 / 2 



AO 



7T 3 frV 3 /3 9 / 2 

where the parameter £ = p,(3 is given implicitly by 
1 



N {i) = 



2hw> 



s /3(«-M) -I- 1 



(25) 



(26) 



with a the chemical potential. 

In Figure 2, the behavior of the dimensionless quantity 
v(N,T) = (V(N,T))^/h7M/(4na N) as a function of 



the temperature is exemplified for N = 10 5 . We can 
observe that although the interaction may be negligible 
for fcgT > 6f its relevance increases dramatically for 



V. DISCUSSION AND CONCLUSIONS. 

We have studied the role of the s— wave contact inter- 
action between Fermi atoms in different states confined 
in a trap by using a simple variational approach, where 
the interaction effects on the variational wave function 
are resumed in an effective frequency. We have found 
closed expressions for the overlap integrals that admit a 
fast convergence and allow a quite direct evaluation of 
pairing energies and total interaction energies. 

Concerning the pairing energies between atoms with 
the same trap energy er = hiu'(n v + 3/2) we determine 
the region of relevance: they become maximum for atoms 
with equal quantum numbers and in states which oscil- 
late in a given axis with excitation number m = n v while 
remain in the minimum oscillating state in the other di- 
rections rij = j 7^ i. For a spherical trap, there is no 
privileged axis of oscillation. Meanwhile, a direct exten- 
sion of the results here obtained, shows that for a com- 
pletely asymmetric trap, oscillations in the direction of 
the highest effective frequency exhibit a higher pairing 
energy. These results can be compared with those found 
for the isotropic trap in spherical coordinates 0] . There 
the individual wave functions have a well defined angular 
momentum I, and it is natural to consider paired state 
functions within a given shell n v that are eigenfunctions 
of the total angular momentum L. In that case the great- 
est pairing energies correspond to n v — I and L = 0. 
These paring energies do not decrease as n v increases. 
For our much simpler two-particle states, 5 max decreases 
as n v . However, our results have a direct translation for 
anisotropic traps, while the relevance of angular momen- 
tum coupled states is highly dependent on the spherical 
'l^ymmetry. Thus, it is important to search for symmetry 
operators that could permit the definition of more tightly 
bound two-particle states for anisotropic traps. 

We have found that the contribution of the contact in- 
teraction potential to the total energy of the fermionic 
ground state behaves as a polynomial, Eq. ^| Being 
valid for both few and many atoms, this expression shows 
the increasing relevance of the interaction terms with 
respect to the trap energies as the number of particles 
grows. In the large N region, the variational equation to 
determine the effective trap frequency relates the varia- 
tional parameter with the geometric factor Rp- Mean- 
while the stability criteria Q restricts the admissible val- 
ues of the coupling constant g. Beyond this limit both 
the model potential and the simple trial wave function 
have not a clear justification. As an extension to finite 
temperature we elaborate on a possible form of the inter- 
action energy expectation value and find that it becomes 
increasingly relevant for fc^T < ep as T diminishes. 
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The approximations made in this work could be valid but 
not just for dilute gases but also for denser systems by 
taking ao as an effective scattering length. A more pre- 
cise study can be made by solving the Hartree-Fock equa- r* _, , f n 

tions of the system using the contact interaction -5] or J dy (cos y+l) m ~ k (cos y - l) k = (-1) 2 m J cos V»- %/2) sin (j 

working with a different model potential and more so- (31) 

phisticated trial wave functions [14| . The computational anc j p^j j 

load of such an approach restricts actual calculations to 

no more than hundreds of particles. Here, we have shown 

explicit results up to 10 6 particles. The possibility of ^ 2 

incorporating correlation effects, that is, working with f sm M-i x cqs^ -1 x dx = —B(— — ) (32) 

more than one Slater determinant for each type of atoms J 2 2 ' 2 

as well as other forms of the effective interaction potential 
is under way. 
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VI. APPENDIX. 



In this appendix, integrals of the form 

/oo 
-oo 



(27) 



Besides llf 



B(r,s) = 



r(r)r(s 



r 



with H n the Hermite polynomials are evaluated. A closed 
expression of these integrals may be obtained employing 
the Rodrigues formula for Hermite polynomials. This 

leads to expressions that contain summations of terms We can thus conclude that 
with alternating signs. Here, we look for an expression 
in terms of finite fast convergent summations. To that 
end, consider the fact [13 



i-^, r(n+l/2) = ^-(2n-l)!! (34) 



Besides 



2 (-1)"(«!) 2 2" f« 



tt (2n)\ 



(cosy) n H 2n (x(l-secy) 1/2 ). 

(28) 



H(n,m) = nlmlJ^ ^ 
" fc=0 



" {r " , " > 4 fc (2m -2k- l)!!(2fc - l)!!(2n -2k- 



(m-k)\(n-k)\(2k)\ 

and the expectation value of the one-dimensional contact 
interaction is then given by 



f 

J — t 



e~ px2 H 2m (bx)H 2n (cx)dx = (2m)\(2n)\J- 



k=0 



(m-k)\(n-k)\(2k%^ 



,b 2 — p^m-k /(? — P 



(6(f = J- 

^ — jjyi-k ^<Z0C^2k V 
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1 



h ^2 n n\ v / 7r2 m m! 



H(n, m) 



so that 
H(n, m) 



/2tt 



Moj x - 

h 



k=0 



(2n- 2k- 1)!! (2m -2k- 1) 
(2n- 2k)\\ (2m- 2k)\\ 



(n\y(m\y fn 



E 

k=0 



(-2) 
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(m-k)\(n-k)\(2k)\ 



dy(cosy + l) n - fe (cosy-l) fc / dy'(cosy + l) m - k (cosy - 1) 



6 



[1] I. F. Silvera and J. T. M. Walraven, J. Appl. Phys 52, 
2304 (1981). 

[2] D. A. Butts and D. S. Rokhsar, Phys. Rev. A 55, 4346 

(1997) . 

[3] M. Houbiers, R. Ferwerda, H.T.C. Stoof, W. I. 

McAlexander, C. A. Sackctt, and R. G. Hulet, Phys. Rev. 

A 56, 4864 (1997). 
[4] J. Schneider and H. Wallis, Phys. Rev. A 57,1253(1998). 
[5] G. M. Bruun and K. Burnett, Phys. Rev. A 58, 2427 

(1998) . 

[6] G. M. Bruun, and C. W. Clark, Phys. Rev. A 61, 061601 

(2000) . 

[7] B. de Marco and D. S. Jin, Science 285, 1703 (1999); A. 

G. Truscott, K. E. Strecker, W. I. Mc Alexander,G. B. 

Patridge, and R. G. Hulet, Science 291, 2570 (2001); F. 

Schreck, L. Khaykovich, and K. L. Corwin, G. Ferrari, 

T. Bourdel, J. Cubizolles, and C. Salomon, Phys. Rev. 

Lett. 87, 080403 (2001). 
[8] S. R. Granade, M. E. Gehm, K. M. O'Hara, and J.E. 

Thomas, Phys. Rev. Lett. 88, 120405 (2002). 
[9] H. Heiselberg, Phys. Rev. A 63, 043606 (2001). 
[10] M. E. Gehm, S. L. Hemmer, S. R. Granade, K. M. 

O'Hara, and J. E. Thomas, Phys. Rev. A 68, 011401 

(2003). 

[11] R. Roth and H. Feldmeier, Phys. Rev. A 64, 043603 

(2001) . 

[12] E. R. I. Abraham, W.I. McAlexander, J. M. Gerton, R. 
G. Hulet, R. Cote, and A. Dalgarno, Phys. Rev. A 55, 
R3299 (1997). 

[13] H. Heiselberg, B. Mottelson, Phys. Rev. Lett. 88, 
190401(2002). 

[14] J. Carlson, S. -Y.Chang, V. R. Pandharipande, and K. E. 
Schmidt, Phys. Rev. Lett. 91, 050401(2003). 

[15] I.S.Gradshteyn and I. M. Ryzhik, Table of Integrals Se- 
ries and Products, Academic Press, New York, 1980. The 
relevant equations are 7.389, 3.621.5, 8.38.4.6 and 8.35. 

[16] A. P. Prudnikov, Yu A. Brychkov, O. I. Marichev, Inte- 
grals and Series, Vol. II, Gordon and Breach, New York 
(1986). The relevant equation is 2.20.16.4. 



0.07- 



0.06- 



5 




50 100 150 



20- 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 



k B T/E F 

FIG. 2: Behavior ol the interaction v(N,T) for N = 10 5 
atoms as a function of the temperature. 



FIG. 1: Behavior of the dimensionless pairing interaction 5 nv 
as a function of the energy level quantum number n v . The 
allowed region is delimited by its maximum 5 ma x (top curve) 
and minimum 8 min (bottom curve). 



